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ABSTRACT 



o. 

Chemistry plays an important role in the structure and evolution of protoplanetary 
disks, with implications for the composition of comets and planets. This is the first 
of a series of papers based on data from DISCS, a Submillimeter Array survey of the 
chemical composition of protoplanetary disks. The six Taurus sources in the program 
(DM Tau, AA Tau, LkCa 15, GM Aur, CQ Tau and MWC 480) range in stellar spectral 
type from Ml to A4 and offer an opportunity to test the effects of stellar luminosity 
on the disk chemistry. The disks were observed in 10 different lines at ~3"resolution 
and an rms of ~ 100 mJy beam -1 at ~ 0.5 km s _1 . The four brightest lines are CO 
2-1, HCO + 3-2, CN 233/4/2 — 1 2 2/3/1 and HCN 3-2 and these are detected toward all 
sources (except for HCN toward CQ Tau). The weaker lines of CN 222 — li 1, DCO + 3— 
2, N2H 4 " 3-2, H2CO 3o 3 — 2o 2 and 4i 4 — 3i 3 are detected toward two to three disks each, 
and DCN 3-2 only toward LkCa 15. CH3OH 42 1 — 3i 2 and 0C3H2 are not detected. 
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There is no obvious difference between the T Tauri and Herbig Ae sources with regard 
to CN and HCN intensities. In contrast, DCO+, DCN, N 2 H+ and H 2 CO are detected 
only toward the T Tauri stars, suggesting that the disks around Herbig Ae stars lack 
cold regions for long enough timescales to allow for efficient deuterium chemistry, CO 
freeze-out, and grain chemistry. 

Subject headings: protoplanetary disks; astrochemistry; stars: formation; ISM: molecules; 
techniques: high angular resolution; radio lines: ISM 



Introduction 



Protoplanetary disks link protostars and (extra-)solar planetary systems physically and chem- 
ically. Understanding the chemical composition and evolution of disks thus provides constraints on 
the nature of molecules incorporated into planetesimals and planets. A variety of simple species, 
including the organic mole cules CN, HCN, and H2CO, have been detected toward s a handful of 



disks in unresolved studies foutrev et al.l 119971 : iThi et al. 



2004 



Kastner et al.ll2008l ). suggestive of 



an active ch emistry. However, apart from CO and to some extent HCO + , the chemistry is poorly 
constrained (jPietu et al.l 120071 ). 



Observations of the earlier stages of star formation and of comets, the possible remnants of our 
own protoplanetary disk, re yeal a chemistry rich in simple and complex organics up to HCOC9H 5 



and (CH20H)2 in size (e.g. Ivan Dishoeck et al. 



1995 



Crovisier et al.l 120041 : iBelloche et all |200J). 



Pre-biotic pathways to chemical complexity thus exist. If theses pathways are active in disks or 
if the observed cometary complexity is instead a fossil remnant from earlier stages remains to be 
shown. Recent experiments suggest that the combination of icy gra in mantles and UV irradiation 
efficiently produces the complex molecules found ar ound protostars (lOberg et al.ll2009l). and surel y 
appropriate conditions are common in disks as well (jvan Zadelhoff et al.ll2003l : iHersant et al.l l2009). 



The chemistry in disks is significantly more difficult to probe than in protostellar cores be- 
cause of their order-of-magnitude smaller angular size, which necessitates the use of (sub-) millimeter 
arrays t o resolve the chemis try of the bulk of the disk material. The first observations of disk chem- 
istry bv lDutrev et al.l (J1997I ) also revealed lower gas-phase abundances of most mole c ules c ompared 
to protostars. From a slightly larger sample consisting of four sources IThi et al.l (J2004I ) showed 
that protoplanetary disks seem to generally contain orders of magnitude lower gas abundances 
compared to protostellar cores. This is reproduced by models of disks with a combination of 
freeze-out onto grains t oward the disk midplane and photodissociation in the disk atmosphere (e.g. 
Aikawa fe: Herbstlll999i ). Gas phase molecules are only expected to be abundant in an intermediate 
zone that is warmer than common ice sublimation temperatures, but sti ll deep enough into the 



disk to be partly p rotected from stellar and interstellar UV irradiation (e.g. lAikawa fc Herbstl ll999: 
Bergin et al.ll2003l ). The low molecular abundances result in weak emission that requires long inte- 
gration times at all existing (sub-) millimeter facilities. The investment required for interferometers 



in large part explai ns the small number of resolved chemistry studies of any species more compli 



cated than HCO+ (jQi et al J 12003 : iDutrevi 12004 : iDutrev et al J 120071 : iQi et aU 12008 : iHenning et al 
20ld ). 



Despite these impediments, single dish studies of the T Tauri stars DM Tau, GG Tau, LkCa 15 
TW Hya, V4046 Sgr and the Herbig Ae stars M WC 480 and HP 163296 have prov i ded some con- 



strain ts on the chemistry of protoplanetary disks (JDutrev et al.lll997l : iThi et alj|2004l : iKastner et al 



20081 ). The species CN, HCN, DCN, HNC, CS, C 2 H, H 2 CO, HCO+ and DCO+ have been detected 
toward at least one of these objects, with CN/HCN ratios that only can be explained by high 
UV or X-ray fluxes penetrating into the disk - CN is a photodissociation product of HCN. The 
variations in molecular abundances among different systems are significant; H2CO is only detected 
toward DM Tau and LkCa 15, DCO + toward TW Hya, and HCN toward all T Tauri stars but 
not toward either of the Herbig Ae stars. Suggested reasons for these variations include higher 
photodissociation rates and a lack of cold chemistry products toward the more luminous Herbig Ae 
stars, different stages of grain growth in different disks, and different disk structures. Differences 
in the chemistry preceding the disk stage may also play a role. 

Resolved s tudies of disk chemi stry are few but intriguing. Using the IRAM Plateau de Bure 
Interferometer, IDutrev et al.1 (J20071 ) reported low signal-to-noise N 2 H + detections t oward DM Tau 



and L kCa 15 and an upper limit toward MWC 480. Within t he same project ISchrever et al. 



(J2008l ) tentatively detected HCN toward AB Aur. More recently IHenning et al.l (|2010l ) resolved 
C 2 H emission toward DM Tau and LkCa 15, while no C 2 H was detected toward the more luminous 
MWC 480. The difference between the T Tauri stars and the Herbig Ae star was expl ained by 



a com bination of high UV and low X-ray fluxes toward MWC 480. Using the SMA, IQi et al. 



(J20081 ) spatially resolved the HCO + and DCO + emission toward TW Hydrae, revealing different 
distributions of these species - DCO + is relatively more abundant at increasing radii out to 90 AU, 
consistent with its origins from cold disk chemistry. HCN and DCN were both detected as well, but 
provide less constraints on the chemistry because of the weak DCN signal. Overall, the combination 
of small samples of diverse sources and even fewer resolved studies have so far prevented any strong 
constraints on the main source of chemical diversity in protoplanetary disks. 

With the DISCS (Disk Imaging Survey of Chemistry with SMA) Submillimeter Array legacy 
program, we aim to produce a resolved, systematic survey of chemistry toward protoplanetary 
disks spanning a range of spectral types and disk parameters. The targeted molecules are the 
simple species that previous studies suggested may be detectable (CO, HCO + , DCO + , CN, HCN, 
DCN, N 2 H + , CsH 2 , H 2 CO and CH3OH). The initial survey contains six well known disks in the 
Taurus molecular clouds (DM Tau, AA Tau, LkCa 15, GM Aur, CQ Tau and MWC 480) with 
central stars that span spectral types Ml to A4. The disk sample is presented in $2] with special 
attention to the properties that may affect the disk chemistry such as stellar luminosity, accretion 
rates, disk size, disk structure and dust settling. The spectral set-ups and observational details are 
described in $3l The channel maps toward one of the richest sources, moment maps for all disks 
in the most abundant species, and spectra of all detected lines toward all sources are shown in §Q 



The detection rates as well as source-to-source variations are discussed in section §5\ followed by 
qualitative discussions on the origins of the observed chemical variations. 



2. The Disk Sample 

2.1. Selection criteria 

The Taurus DISCS sample of protoplanetary disks was chosen to assess the impact of spectral 
type or stellar irradiation field on the chemistry in the disk. The target systems span the full range 
of stellar luminosities among the T Tauri and Herbig Ae stars associated with the Taurus molecular 
cloud. Table Q] lists the stellar properties. The stellar masses in the sample range between 0.5 and 
1.8-2.3 Mq, corresponding to luminosities that span almost two orders of magnitude. If stellar 
luminosity is the main driver of the outer disk chemical evolution, then this sample should display 
a range of chemical behaviors. As discussed below there are other sources of radiation that may 
affect the chemistry as well, especially accretion luminosity and X-rays. 

The sample is biased toward disks of large angular extent, since disks smaller than a few 
hundred AU are not spatially resolved by the SMA in the compact configuration. The observations 
are not sensitive to gas inside of 100 AU, which entails that there is a large gap in radii between 
these millimeter observations and infrared disk chemistry observations that typically probe the 
inner disk out to a few AUs. The sources were selected from disks previously mapped in CO and 
thus may be biased toward gas-rich disks. Furthermore only disks clearly isolated from the parent 
cloud emission are included to reduce confusion and ensure that the detected molecules reside in 
the disks. Known harborers of organic molecules were favored (DM Tau, LkCa 15 and MWC 480), 
but the sample also contains disks that have only upper limits or that have not been investigated 
for molecular lines other than CO. The focus on a single star forming region allows us to probe 
disks of similar ages and likely similar chemical starting points, reducing the sources of chemical 
variation. 



2.2. Star and Disk Properties affecting Disk Chemistry 

Stellar luminosity, accretion luminosity, X-rays, the interstellar irradiation field, disk geometry, 
disk gaps and holes, dust settling and dust growth may all affect the chemistry in the disk. The 
sample characteristics in terms of these properties are discussed in this section. 

Observations of high CN and HCN abundances toward protoplanetary disks reveal that a chem- 
istry controlled b y far-ultraviolet (FUV) or X-ray radiation or both must contribute to the observed 



abundances (e.g. Ivan Zadelhoff et alJbOQll : iThi et all 12004 ). FUV radiation below 2000 A affects 



the chemistry by directly heating the gas in the disk surface layers, in limiting molecular abun- 
dances via photodissociation, increasing the amount of photochemistry products such as CN and 



liberating frozen species via photodesorption. The nature of the dominating source of radiation is 
however unknown. The quiescent stellar luminosities in the sample range from 0.25 to 11.5 L Q , 
and if quiescent FUV radiation controls the photodissociation rate in the disk surface, then there 
should be a clear trend in the CN/HCN ratio between the low-luminosity T Tauri stars and the 
order of magnitude more luminous Herbig Ae stars. However, Kurucz (1993) stellar atmosphere 
models show evidence that stars with a spectral type later than F do not have significant stellar 
continuum <2000 A above that generated by accretion, and that it is only for A stars that stellar 
UV becomes more significant than accretion luminosity for the FUV field. This is confirmed by 
FUV observations tow ard the T Tauri star TW Hya, which is dominated by line emission generated 
from accretion shocks (JHerczeg et al.ll2002l ). 



The accretion FUV spectrum is dominated by line emission and especially Ly-ce e mission, which 



resul ts in preferential HCN dissociation and therefore boosts the CN/HCN ratio (JBergin et al. 



2003) beyond what is expected from a UV chemistry dominated by continuu m radiation. Towar d 
T Tauri stars, the FUV flux is expected to scale with the accretion rate (e.g. ICalvet et al.l 12004 ) . 
The accretion rates proba bly vary over time however, as has been observed toward TW Hya, where 
Alencar &; Batalhal ( 2002 ) found mass accretion rates between 10 -9 - 10 -8 M yr _1 during a 
year and smaller variations on timescales of days. The measured accretion rate variability among 
the T Tauri stars in this sample are all within this range and it is unclear whether the average 
accretion luminosity v ary significantly am ong the sources. In general more massive systems have 
higher accretion rates (jCalvet et al.ll2004l ) and it can therefore be expected that the FUV flux from 
accretion will be higher toward the intermediate mass Herbig Ae stars in the sample. As mentioned 
above, for the early type stars (e.g. A) the stellar continuum also adds a significant contribution to 
the FUV radiation field. Thus, even if accretion rather than quiescent luminosity controls the UV 
chemistry, one would expect to observe a higher CN/HCN ratio toward the more massive stars. 

Observations of the FUV field in T Tauri systems find that accretion produces UV fluxes 
that are a f ew hundred to a th ousand times stronger than the Interstellar Radiation Field (ISRF) 
at 100 AU (JBergin et al.1 12004 ) . The ISRF may still be important at large radii, however. The 
external irradiation field is presumed to be constant toward the Taurus sample of sources, i.e. 
none of the disks are close to any O or B stars, but its impact may differ between the different 
classes of sources; the ISRF may play a larger role in driving the chemistry for the low luminosity 
objects, thus acting as a chemical equalizer. A fourth possible driver of the disk chemistry is X 



rays, which is predicted to be important for the ionization fraction in the disk (JMarkwick et al 



2002). X-rays are mainly atte nuated by gas, while continuum UV photons are quickly absorbed 
by dust (JGlassgold et al.lll9971 ). X-rays can therefore penetrate deeper into the disk compared to 
UV rays, and may be a main driver of both ion chemistry and molecular dissociation. This may 
result in observable differences in the molecular distribution if the chemistry is driven by X-rays 
instead of UV radiation. For individual objects X-ray measurements are notoriously variable and 
assessing their impact observationally may be possi ble only through mon itoring; of the X-ray flux 
and the chemical variation toward a single system (JGlassgold et al.ll2005l ). On average the X-ray 



fluxes seem higher toward T Tauri stars compared to Herbig Ae stars and t hus ion-driven chemical 
reactions may be faster in disks around T Tauri stars (JGiidel fc Nazel 120091 ) . 



Regardless of whether the quiescent luminosity of the central star is the main source of energetic 
radiation, it may still control the disk temperature and thus the temperature sensitive chemistry, 
e.g. deuterium fractionation efficiency probed by the DCO + /HCO + ratio. The stellar continuum 
photons are expected to be the primary agent for heating grains in the disk, particularly toward the 
midplane, setting the overall reservoir of warm or cold grains in the outer disk to which the SMA is 
most sensitive. Both the level of gas phase depletion and the chemistry dependent on CO depletion 
may then be regulated by the stellar continuum flux. For example a tracer of CO fr eeze-out such 



as th e N2H + /HCO + ratio is expected to be higher toward low-luminosity systems (JBergin et al 
20021 ). 



The disk structure and dust properties determine how much of the stellar radiation is inter- 
cepted by the disk and the penetration depth of that radiation. These disk properties may be 
equally important to the the strength of the radiation field for the chemical evolution in the disk. 
Table [2] list the sample disk characteristics from previous CO and continuum observations and 
modeling. Position angles and inclinations do not affect the chemistry intrinsically, but the disk in- 
clination affects which parts of the disk are observed, and thus our view of the chemistry. The sizes 
of the disks are described in terms of CO gas disk radii and range from 200 to 890 AU. Disk masses 
are less well constrained, since they are derived from dust emission assuming a dust-to-gas ratio. In 
most studies the canonical interstellar dust-to-gas ratios of 1 to 100 is used. The actual dust-to-gas 
ratio may be different because of dust coagulation and photoevaporation, and also variable among 
the sources. Still it seems clear that two of the disks, AA Tau and CQ Tau, are significantly less 
massive than LkCa 15, GM Aur and MWC 480, while DM Tau falls in between (Table [2]). This 
difference in disk mass may affect their relative abundances of different species, since more massive 
disks are expected to contain more cold material per unit of incoming irradiation. 

Three of the sample disks are so-called t ransition and pre- t ransition disks (DM Tau, GM Aur 



and L kCa 15) with large inner holes or gaps (jCalvet et al.ll2004l : lEspaillat et al.1 120071 : iDutrey et al. 



2008). From a Spitzer survey of disk chemistry ma ny transition disks, including GM Aur and LkCa 
15 have CO gas in the disk hole (jSalvk et al.ll2009i ). They do however lack emission from HCN and 



C2H0 transitions in the inner disk that are strong toward classical T Tauri stars (jPascucci et al. 



2009). It is unclear whether this is a chemistry or gas-mass effect. It is also unknown whether the 
chemistry in the outer regions of these transition disks differ significantly from classical T Tauri 
disks, although it is curious that these systems are the most chemically rich in terms of the outer 
disk seen to date. This may be explained by increased radiation fluxes; more stellar radiation may 
reach the outer disk the larger the hole. In addition, large holes may be a tracer of overall grain 
growth and thus of increased UV penetration depth in both the inner and outer disk. 

Grain properties can be traced by the millimeter slope of the Spectral Energy Distributions 
(SEDs), parameterized by the power-law index of the opacity spectrum, /3, and this is the last disk 



characteristic listed in Table [21 The index j3 is predicted to decrease with grain growth and the 
disks in the survey all have /3 below the expected 1.7 for interstellar grains ([Andrews fc Williams 
20071 ). The (3 estimates are however different between different studies, which makes it difficult to 
conclude on an order of dust coagulation among the disks, but A A Tau seems to have a significantly 
lower /3 compared to the other T Tauri disks. Most of the sources h ave also been obse rved by 
Spitzer in studies that constrain the grain properties in the inner disk (JFurlan et alJl2009l ). These 
measurements do not provide any straightforward constraints on the grain properties in the outer 
disk probed by the SMA, however. 

Finally the disk structure, e.g. the amount of flaring, affects the amount of stellar light inter- 
cepted by the disk. The disk structure, parameterized by the dust scale heigh t, can be constrained 



by modeling the SED, though degeneracy is a problem ([Chiang et al. 



2001 



As the dust settles 

towards the midplane the dust scale he ight decreases comp ared to the gas scale height. Five of 
the sample disk SEDs were modeled by IChiang et al.l (J200ll ) , who concluded that MWC 480 and 
LkCa 15 are more settled than CQ Tau, GM Aur and AA Tau, i.e. in the disks of MWC 480 and 
LkCa 15 the upper disk layers are significantly depleted in micron-sized dust grains. This has been 
confirmed in more recent stud i es that find signific antly less settling toward GM Aur compared to 



LkCa 15 (JEspaillat et al. 



2007 



Hughes et al.1 12009 ). 



Observations 



3.1. Spectral Setups 

Two frequency setups per source were selected to cover 4 to 8 spectral lines in each setup, 
at the same time providing continuum observations. The targeted molecules are DCO + and DCN 
(probing deuterium chemistry), CN and HCN (probing photochemistry), HCO + and N2H + (ions), 
H2CO and CH3OH (potential grain chemistry products), 0C3H2 (carbon-chain chemistry) and CO 
(the disk kinematic tracer). Tables [3] and H] summarizes the two spectral setups, centered at 1.1 
mm and 1.4 mm, respectively. 

The SMA correlator covers ~ 4 GHz bandwidth in each of the two sidebands using two Inter- 
mediate Frequency (IF) bands with widths of 1.968 GHz. The first IF band is centered at 5 GHz, 
and the second IF band is centered at 7 GHz. Each band is divided into 24 slightly overlapping 
"chunks" of 104 MHz width, which can have different spectral resolution. For each spectral setting, 
the correlator was configured to increase the spectral resolution on the key species with 128-256 
channels per chunk, with the exception of one H2CO line observed with 64 channels. Chunks 
containing weaker lines were then binned to obtain higher signal to noise, while still recovering 
sufficient kinetic information. The remaining chunks were covered by 64 channels each and used 
to measure the continuum. The continuum visibilities in each sideband and each IF band were 
generated by averaging the central 82 MHz in all line- free chunks. 



3.2. Data Acquisition and Calibration 

The six disks were observed from 2009 November through December with the compact config- 
uration of the Submillimeter Array (SMA) interferometer (Ho et al. 2004) at Mauna Kea, Hawaii. 
LkCa 15 was also observed on 2010 March 22nd with the subcompact configuration for more short- 
spacing data to improve the signal-to-noise of deuterated species detection. For each observation, 
at least six of the eight 6 m SMA antennas were available, spanning baselines of 16-77 m (Table 
[5]). The observing sequence interleaved disk targets and two quasars in an alternating pattern. De- 
pending on their proximity to the disk targets and fluxes at the time of the observations, a group 
of three quasars was used as gain calibrators: 3C 111, J0510+180 and 3C 120. The observing 
conditions were generally very good, with T225GHZ ~ 0.04 — 0.1 and stable atmospheric phase. 

The data were edited and calibrated with the IDL-based MIR software packagqj. The bandpass 
response was calibrated with observations of Uranus and the bright quasars available (3C 454.3, 3C 
273, 3C 84). Observations of Uranus provided the absolute scale for the calibration of flux densities 
for all compact tracks and Vesta for the one LkCa 15 sub-compact track. The typical systematic 
uncertainty in the absolute flux scale is ~10%. Continuum and spectral line images were generated 
and CLEANed using MIRIAD. MIRIAD was also used to calculate synthesized beam sizes and 
position angles, which are listed in Table [6] for each setting and source. 



4. Results and simple analysis 

The first section presents the extracted spectra for all molecules detected toward at least one 
disk, channel maps toward DM Tau, which contain some of the strongest detections of the weaker 
lines in the survey and finally disk images and moment maps for strong molecular lines toward all 
sources ( §4.1|> . The subsequent sections provide more details on the observed CN/HCN ratios, ions 
and deuterated molecules and their variation among the surveyed disks. 

In general, we do not attempt to estimate column densities of molecules, but rather present the 
data in terms of integrated fluxes and flux ratios. The line fluxes and column densities are of course 
related, but the fluxes also depend on excitation temperatures and opacities. Determining molecular 
column densities therefore requires knowledge of the disk structure and the spatial distribution of 
molecules. Even for optically thin lines, estimated column densities may be off by an order of 
magnitude or more if the wrong emission regions/temperatures are adopted. For optically thick 
lines the estimates become even more uncertain and single dish studies including 13 CO and H 13 CQ + 



have shown that the emission from 12 CO and H 12 CO + is mostly optically thick ( Dartois et al. 



2003 



Thi et alj|2004l ). A proper derivation of column densities therefore requires detailed chemical 



modeling of the disks, which we anticipate for a future DISCS publication. An exception from only 



1 http://www.cfa.harvard.edu/~cqi/mircook.html 



reporting fluxes is made for the deuterium fractionation, where it is useful to derive abundance 
ratio limits assuming the same emission regions of molecules and their deuterated equivalents, to 
compare with previous studies. We also estimate the HCO+/CO ratio toward CQ Tau for a rough 
consistency check with previous observations. 



4.1. Overview of Taurus disk chemistry 

Figure [T] shows the extracted spectra of the targeted molecules toward DM Tau, AA Tau, LkCa 
15, GM Aur, CQ Tau and MWC 480. The spectra are also tabulated in Table The four spectral 
lines, CO 2-1, HCO + 3-2, HCN 3-2 and CN 3-2, expected to be strongest from previous studies 
were detected in all six disks except for HCN toward CQ Tau, with an order of magnitude variation 
in integrated fluxes between the different disks. Spectral lines from N2H + and H2CO were detected 
toward three disks, DM Tau, LkCa 15 and GM Aur, with tentative detections toward AA Tau. 
Lines from the two deuterated molecules, DCO + and DCN, were detected toward LkCa 15 and the 
DCO + 3-2 line was also detected toward DM Tau. The spectral lines from CH3OH and C-C3H2 
were not detected toward any of the disks. 

Figure [2] shows the velocity channel maps for all detected spectral lines toward DM Tau. The 
peak intensities vary between ~9 Jy beam -1 for CO and ~0.1 Jy beam -1 for the weaker H2CO 
line. All species, except for H2CO, are apparent in multiple velocity channel maps, and the stronger 
detections clearly follow the velocity pattern indicative of a disk in Keplerian rotation. The relative 
intensity of different species in different velocity channels varies significantly, indicative of differences 
in emission regions between different lines. 

Fig. [3] shows integrated intensity and first moment maps derived from the channel maps for 
CO, HCO + , HCN and CN for each disk. In addition, the first column in Figure [3] shows millimeter 
continuum SEDs compiled from the literature □ together with new measurements at 218 GHz and 
267 GHz (Table [8]), which show good agreement for all of the sources. The dust continuum flux 
densities at ~267 (and 218) GHz vary by about a factor of four among these sources, with MWC 
480 the strongest at 430 mJy and AA Tau the weakest at 110 mJy. Note that there is not a 
one-to-one correspondence between the strength of the 218 or 267 GHz dust continuum emission 
and the CO 2-1 line emission. 

The first-moment maps show the Keplerian rotation pattern of the disks, most clearly in the 
CO 2-1 and HCO + 3-2 emission. The rotational velocity pattern is also present to some extent in 
HCN and CN emission, which suffer from lower signal-to-noise and, for CN, the blending of the 
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spectroscopic triplet. 

Line spectra are extracted from the channel maps using elliptical masks produced by fitting 
a Gaussian profile to the CO integrated intensity maps toward each source to obtain major and 
minor axes and positions angles (listed in Fig. [3]). The size of the mask is scaled for each line, to 
optimize the signal-to-noise without losing any significant emission, such that the major and minor 
axes are between 2-<r and 1-a of the Gaussian fitted to the CO emission - a 2-a radius corresponds 
to ~1.7 x the radius at full-width-half-maximum (FWHM). All masks are large enough to not 
cut out any emission, i.e. the chosen masks results in no significant decrease in integrated line 
intensity compared to integrating over the full CO disk. Applying these masks ensures that only 
disk emission and minimal noise are included in the spectra. The derived CO full-width-half- 
maxima agree reasonably well with the disk sizes from the literature compiled in Table [21 i.e. all 
semi-major axes are within 1" of the previously observed or derived CO radii, where resolved data 
exist. In addition the position angles agree within 10° except for the barely resolved disk of CQ 
Tau. The integrated spectra from CQ Tau are however not significantly affected by the choice of 
mask shape. 

The resulting spectra in Fig. [1] are used to derive the total intensities listed in Table The 
2-0" upper limits are calculated from the rms when the spectra are binned to a spectral resolution 
of 3.3-4.3 km s , the minimum resolution for resolving any lines, multiplied by the full width 
half maximum of the CO 2-1 line toward each source. This approximation is supported by the 
similarity in the line widths for different transitions toward the same disk in Fig. [H Because of 
variable observing conditions and integration times, the flux upper limits vary between 0.09 and 
0.78 Jy km s _1 per beam. Most upper limits are however lower compared to detected fluxes toward 
other sources. Specifically the non-detections toward MWC 480 of DCO + , N2H + and one of the 
H2CO lines are up to a factor of two lower than the detected fluxes toward the T Tauri stars DM 
Tau, GM Aur and LkCa 15. 



4.2. CN/HCN flux ratios 

As discussed above, deriving molecular abundances from disk spectra is complicated by the 
uncertainties in disk structures and the accompanying uncertainties in the emission conditions of 
the different molecules. Line flux ratios can however be used as a proxy for abundance ratios when 
the emission is optically thin and the upper energy levels are similar, canceling out the temperature 
effect, if the same emission region for the two molecule s can be assumed. M odels suggest that the 



emission regions are different for CN and HCN, however (jJonkheid et al.ll2007l ). It is still informative 
to compare the flux ratios between different sources if it can be assumed that the relative emission 
regions of the CN and HCN is stable between different disks, i.e. that CN is always present in a 
warmer layer closer to the surface compared to HCN. Then changes in the CN/HCN flux ratio can 
still be used to trace changes in abundance ratios if the emission is optically thin. 
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Thi et al.l (|2004l ) suggested that HCN and CN line emission from protoplanetary disks may 
be somewhat optically thick from comparisons of HCN and H 13 CN line intensities toward LkCa 
15. This survey does not currently include any rare isotopologues, but the CN optical depth can 
be estimated from the relative integrated intensities of the 23-I2 triplet transition and the ~10 times 



wea ker 2^-1 1 singlet tra nsition using the line strengths from CDMS (http://www.astro.uni-koeln.de/cdms/catalog 



and lMiiller et al.ll200ll ) at any of the reported temperatures between 9 and 300 K (the transitions 
have the same excitation energy). Toward DM Tau and AA Tau the observed CN line ratios are 
consistent with optically thin emission. Toward LkCa 15 and MWC 480, the emission from the 
CN triplet underestimates the CN abundance by factors of 1.3 and 1.6, respectively, indicative of 
somewhat optically thick CN triplet emission. 

Despite the complications introduced by modest optical depth, changes in CN 23-I2 / HCN 
3-2 ratios larger than a factor of 2, assuming similar levels of optical depth for HCN emission, are 
expected to trace real variations in the chemistry as discussed above. The absolute CN and HCN 
integrated intensities range from 0.2 Jy km s _1 toward CQ Tau to 5.5 Jy km s _1 toward LkCa 15. 
Figure H] shows that the variation in CN/HCN line intensities is smaller - all sources have ratios of 
0.8-1.6 except for AA Tau, which has a CN/HCN emission ratio of 2.9. AA Tau is then the only 
significant outlier in terms of CN/HCN flux ratios. 



4.3. Ions: HCO+ (DCO+) and N 2 H+ 

HCO + is often used as a tracer of gas ionization and thus of high energy radiation in disks. The 
high optical depth of HCO + and the lack of rare isotopologues of CO prevent such an analysis at 
prese nt, except for toward CQ Tau where the CO 2-1 emission has been estimated to be optically 



thin (jChapillon et al.ll2008l ). The ratio of column densities for species X and Y, for optically thin, 



resolved emission, can be calculated from 

Ny fT&dv Ql t (T) e W« ^SxMx' U 

where N is the column density, J T^dv is the integrated line emission in K km s , Qrot(T) 
the temperature dependent partition function, E u the energy of the upper level in K, T ex the 



excit ation temperature in K and Sy^i 2 are the line strength and dipole moment (e.g. iThi et al 



2004! ) . Toward the same source, the integrated line flux in Jy and line intensity in K are related by 
T m b[K] oc .F[Jy] x A 2 [mm 2 ]. Using partition functions, level energies and line strengths from CDMS 
and assuming the same excitation conditions for HCO + 3-2 and CO 2-1 the [HCO + ]/[CO] ratio is 
1.0 — 1.3 x 10 -4 toward CQ Tau for excitation temperatures of 18-75 K. 

The N2H" 1 " observations toward DM Tau, LkCa 15 and GM Aur (and tentatively toward AA 
Tau) pro yide unambiguous d etections of this species in protoplanetary disks, confirming previous 



claims bv lDutrev et al.l (|2007l ). Since N2H + and DCO + both potentially trace the chemistry further 
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toward the midplane compared to the more abundant molecules, their ratio provides important 
constraints on the cold chemistry in disks. The DCO + /N2H + line intensity ratio ranges from <0.1 
toward GM Aur to 0.8 toward DM Tau (Fig. \S§. In contrast, there is no significant variation in 
the H2CO/N2H+ ratio between the sources. Assuming that these molecules always reside in the 
colder regions of the disks, these differences in flux ratios suggest relative abundance variations of 
N2H + and DCO + of an order of magnitude between the different sources. 



4.4. Deuteration: DCO+/HCO+ and DCN/HCN 

Because of the high optical depth of the HCO + line emission, the DCO + /HCO + line intensity 
ratio can only be used to derive upper limits on the average HCO + deuteration fraction in the 
disk. T he analysis is f urther complicated by evidence of different emission regions of DCO + and 



HCO + ()Qi et al.ll2008l ). Assuming, however, the same emission region of DCO + and HCO + and 
optically thin emission for both ions, the upper limits on the deuteration fraction is calculated using 
Eq. Q] to vary between 0.32 toward DM Tau, 0.18 toward LkCa 15 and <0.07 toward GM Aur. 
An excitation temperature of 19 K is assumed, but the ratios are only marginally affected by the 
excitation temperature between 10 and 50 K. The variable ratios hint at differences in deuteration 
fractionation between different sources, especially since GM Aur has both lowest DCO + /HCO + 
and DCO + /N2H + ratios, but radiative transfer modeling is needed to confirm this result. 

DCN is only detected toward LkCa 15. The detection is at the >5-cr level, and the moment 
maps shows an almost perfectly aligned velocity field compared to CO and HCO + and the detection 
appears secure. From the CN analysis above and previous single-dish observations, HCN line 
emission is expected to be much less optically thick than HCO + emission and the DCN/HCN ratio 
should provide stricter limits on the deuteration level in the disk. Assuming optically thin emission, 
the same emission region for DCN and HCN and an excitation temperature of ~40 K, the upper 
limit on average deuteration in the disk around LkCa 15 is 0.06, a factor of three lower than the 
estimate from DCO+/HCO+. 



5. Discussion 

5.1. Detection rates and comparison with previous studies 

The reported images and spectra were acquired with only 3-7 hours on source integration, 
with the shorter time spent in the 1.1 mm setting. The detection rate of N2H + and H2CO in the 
1.1 mm setting is then quite remarkable and can be attributed to the advantages of targeting higher 
J lines when probing disks. For comparison, N2H + was previously detected toward DM Tau and 
LkCa 15 through it s 1-0 line using the Plateau de Bure Interferometer, with peak fluxes <0.02 Jy 



(JDutrey et al.1 120071 ). an order of magnitude or more lower than the 3 — 2 line peak fluxes reported 
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here. The 1-0 line in LkCa 15 was also detected bv lQi et al.l (|2003l ) using the Owens Valley Radio 
Observatory Millimeter Array with integrated intensity four times larger than that with the PdBI. 

The agreement with previous single-dish observations is generally good. All targeted molecular 
lines that have been previously report ed in the single dish studies of DM Tau, LkCa 15 and MWC 
480 are also detected with the SMA (|Dutrev et al.lll997l : VThi et alJbooJ : iGuilloteau et alJboOfih . 
Where the same lines have been studied, most integrated in tensities agree . HCN J = 3 — 2 toward 
DM Tau is an exception, where the reported upper limit in iDutrev et al.l (j 19971 ) is a factor of two 
below the intensity observed with the SMA. 



Using their detections, butrev et al.1 (J2007I ) derived [N 2 H+]/[HCO+] ratios of 0.02-0.03 for 
DM Tau and LkCa 15 by fitting the line emission to disk models. Without such modeling we 
can only derive upper limits on the [N2H + ]/[HCO + ] of 0.13-0.19 for the two disks because of 
the HCO + line optical depth. Considering that the HCO + abundance may be underestimated by 
up to a n order of magnitude, the two data sets are consistent. Within the same observational 



program 



Chapillon et al.l (|2008l ) searched for CO 2-1 and HCO + 1-0 emission toward CQ Tau and 



used the data to derive an upper limit on the [HCO + ]/[CO] abundance ratio. Assuming the same 
CO and HCO + distribution and excitation conditions and optically thin CO emission they find 
[CO]/[HCO + ]>4xl0 3 . This is consistent with the abundance ratio of 10 report ed above, which 



i s calc ulated making similar assumptions, but without the detailed modeling in IChapillon et al. 



CQ Tau is by far the most chemically poor of the investigated disks. It is interesting that 
despite the low abundances, the chemistry appears 'normal', the ratios of the integrated intensities 
toward MWC 480 and CQ Tau are the same within a factor of two, including the CN/HCN emission 
ratio. The only difference is that overall the gas toward CQ Tau is probably richer in CN and HCN 
with respect to CO, taking into account the large optical depth of the CO emission toward MWC 
480 (jThi et al.ll2004l ). as might be expected for a smaller disk, completely exposed to UV radiation. 



The upper limits on the [DCO + ]/[HCO + ] abundance ratio of <0. 07-0. 32 found toward the T 
Tauri systems in DISCS are consistent w ith the ratio of 0.035-0.05 observed toward TW Hydrae 
van Dishoeck et al.l[20oi : loi et al.l 120081 ). The better constrained [DCN]/[HCN] ratio of 0.06 to- 



ward LkCa 15 is also consistent with the value of 0.02-0.05 toward TW Hydrae. While this ratio 
may be overestimated by a factor of a few, high levels of deuterium fractionation seems common 
toward T Tauri systems. 

As found in single dish studies, the CN/HCN line intensity ratios toward the Taurus disks 
are high compared to interstellar clouds and cores. The line ratios all fall within the range of 
measurements towa rd other disks, where the total integrated flux ratio of CN/HCN is ~l-5 (see 
Kastner et al.ll2008l . for a compilation). A more quantitative comparison with previous observations 
is difficult without detailed modeling because different studies observed different transitions of CN 
and HCN. 
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5.2. T Tauri vs. Herbig Ae stars 

In agreement with previous studies, we find that the disks surrounding T Tauri stars are more 



chemically rich in spe c ies with strong mm-transitions com 



(jChapillon et al. 



2008 



Schrever et al. 



2008 



Henning et al. 



pared to disks around Herbig Ae stars 



2010 ). The observed chemical poverty 



in the outer disks of Herbig Ae stars has been attributed to the more intense UV field around 
Herbig Ae stars compared to T Tauri stars, which may efficiently photodissociate most targeted 
molecules. In this sample, the most obvious difference between T Tauri and Herbig Ae stars is the 
lack of the cold chemistry tracers N 2 H + , DCO+, DCN and H 2 CO toward CQ Tau and MWC 480, 
while they are detected in 3/4, 2/4, 1/4 and 3/4 of the T Tauri systems. The upper limits toward 
CQ Tau are less informative because of its weak CO emission and observations toward more Herbig 
Ae stars are required to confirm that this difference between disks around low and medium mass 
stars is general. 

In contrast the CN and HCN emission is similar toward the lowest and highest luminosity 
stars, DM Tau and MWC 480, in the sample. CN and HCN emission are modeled to originate 
mainly from the outer layers of the disk (IWillacyll2007l ) and this chemistry thus seems equally active 
toward low- and intermediate-mass pre-main sequence stars. 



5.3. CN and HCN 



CN is a photodissociation product of HCN and the CN/HCN ratio has been put forward to 
trace several different aspect s of the UV field. The C N/HCN ratio is proposed to increase with 
the strength of the UV field (Ivan Zadelhoff et al.ll2003l ). and it will be further enhanced if the UV 
radiation is domi nated by line emiss ion from accretion, since HCN is dissociated by Ly-a photons 
while CN is not (IBergin et al.ll2003l ). Dust settling or coagulation allows radiation to penetrat e 
deeper into the disk, which is also predicted to enhance the CN/HCN ratio (jJonkheid et al.1 120071 ). 



The quiescent UV luminosity increases with stellar mass. There is however no visible trend in 
the emission ratio of CN/HCN with spectral type. In fact, all CN/HCN ratios are the same within 
a factor of two, except toward AA Tau, which has a factor of a few higher intensity ratio. This 
suggests that the CN/HCN ratio is not set by the stellar luminosity though there are complications 
in comparing CN/HCN ratios toward disks around low- and intermedia te-mass stars b ecause of 
potentially different excitation conditions for HCN in the two sets of disks (jThi et al.ll2004l ) . Within 
this sample the CN/HCN ratio also does not trace accretion luminosity; AA Tau has a comparable 
accretion rate to LkCa 15 and among the sources with comparable CN/HCN ratios the accretion 
rate varies by an order of magnitude. AA Tau is reported to have a lower power-law index of the 
opacity spectrum, /?, compared to the other disks, indicative of dust growth and it may be the dust 
properties rather than the stellar or accretion luminosities govern the importance of photochemistry 
in disks. A larger sample that spans a wider range of accretion rates and dust properties is clearly 
required to give a more definitive answer. 
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An additional complication is that the high CN/HCN ratio to ward AA Tau may b e a geometric 
effect. Compared to the other disks, A A Tau is almost edge-on (jMenard et al.ll2003l ). which may 
result in preferential probing of the disk atmosphere compared to less inclined disks. Disk chemistry 
models (Fogel et al. submitted to ApJ) show that CN mainly emits from the disk surface, while 
HCN emission originates further into the molecular disk layer and the more inclined disk may offer 
a viewing angle that is biased toward CN emission. To estimate the effect of disk inclination on 
CN/HCN flux variations then requires a combination of chemical modeling and radiative transfer 
models. 

In terms of absolute flux intensities, the weak CN and HCN emission toward GM Aur compared 
to LkCa 15 and DM Tau stands out. The difference between GM Aur and LkCa 15 may be due to 
the higher accretion rate and intenser FUV field toward LkCa 15. The difference between DM Tau 
and GM Aur is however difficult to explain in terms of UV flux, since DM Tau is a weaker accretor 
than GM Aur. There is so me evidence for significantly more dust settling toward LkCa 15 and DM 



Tau compared to GM Aur (|Chiang et al. 



2001 



Espaillat et al. 



2007 



Hughes et al.ll2009l ). This may 



expose more of the gas in the LkCa 15 and DM Tau disks to high-energy radiation, enhancing the 
photoproduction of CN and HCN as well as the ion chemistry deeper in toward the disk midplane. 



5.4. Cold chemistry tracers 

Lower abundances of DCO + , DCN, N2H + and H2CO toward more luminous stars are qual- 
itatively consistent with our current chemical understanding. DCO + fo rms efficiently from gas 
phase reactio ns with H2D + , which is only enhanced at low temperatures (jRoberts &: Millarl l2000l : 
Willacy! 120071 ) and should be enhanced toward colder disks. Among the T Tauri stars the brightest 
DCO + emission is observed toward the disk around the least luminous star, DM Tau, consistent 
with a higher degree of deuterium fractionation around colder stars. The difference in DCO + line 
flux around GM Aur and LkCa 15 is more difficult to explain. GM Aur has a more massive dust 
disk than LkCa 15 and the two stars have similar luminosities. Naively GM Aur should then be 
surrounded by at least as much cold disk material as LkCa 15. Instead, the upper limit on the 
DCO + flux is a factor of three lower toward GM Aur compared to LkCa 15. There is thus no 
one-to-one correlation between the ratio of disk dust mass over quiescent stellar luminosity and 
DCO + column densities. 

N2H" 1 " forms f rom protonation of N2 by Hjj~ and is mainly destroyed by reactions with CO 
(JBergin et al.ll2002l ). Abundant N2H + is therefore only expected where CO is depleted onto grains 
toward the disk midplane. N2 freezes onto grains a few degrees below CO (jOberg et al.ll2005l ) and in 
cold disks the N2H + abundance should peak in a narrow region where the temperature is between 
the N2 sublimation temperature of ~16 K and the CO sublimation temperature of ~19 K. As long 
as a cold region exists in the disk, the N2H + abundances may be quite independent of the total 
amount of cold disk material. The observations are consistent with this disk abundance structure; 
within the T Tauri sample the N2H + emission only varies by a factor of two, increasing slightly 
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with increasing disk mass. 

The variation in DCO + /N2H + flux ratios over the sample suggest that while both molecules 
trace a cold chemistry, their dependences on the physical environment is considerably different. 
The >8 times higher DCO + /N2H + flux ratios toward DM Tau and LkCa 15 compared to GM 
Aur may be related to the 5-6 times higher fluxes of CN and HCN toward DM Tau and LkCa 15 
compared to GM Aur. This would suggest that both ratios depend on the amount of dust settling 
and that DCO + trace a cold radiation driven chemistry. Considering the ions involved in forming 
DCO + and DCN, it seems reasonable that their formation will be enhanced in regions that are 
irradiated, but not heated by FUV photons or X-rays. To test this hypothesis requires H 13 CO + 
abundances toward both systems (to measure whether the DCO + /HCO + abundance ratio varies 
as well) in combination with a model that simultaneously treat deuterium chemistry and UV and 
X-ray radiative transfer. 

H2CO can form both th rough gas and grain su rface processes. The gas phase process starts 



with CH^ reacting with H2 ( Roberts Sz Millar! 120071 ) and is expected to be at least as efficient in 



disks around low and intermediate mass stars. In contrast, H2CO formation on grains requires 
the freeze-out of CO, which is only efficient at low temperatures. The absence of H2CO toward 
the more luminous stars suggests that the grain surface formation mechanism dominates in disks 
and it is also another indication of the lack of a large cold chemistry reservoir toward disks around 
intermediate mass stars. It also suggests that the organic molecules formed in the protostellar 
stage, where H2CO is common, do not survive in the gas phase in mature disks. 

In summary, all potential tracers of cold chemistry imply the same lack of cold disk material 
aroun d Herbig Ae stars, which is i n agreement w ith a recent survey of CO gas toward Herbig Ae/Be 



stars ([Panic &; Hogerheiidd 12009). In contrast, iPietu et al.l (J2007I ) find that the disk around the 



Herbig Ae star MWC 480 contains large amounts of cold CO gas, below 17 K, indicative of cold 
material outside of 200 AU. At these temperatures CO should not be in the gas phase at all, however, 
since it is below the sublimation point of CO ice. Its presence is a sign of either effi cient mixing in the 



disk or efficient non -thermal ice evaporation, perhaps through photodesorption (jOberg et al. 



2007 



Hersant et al.ll2009l ) . Mixing may drag up material from the midplane on shorter timescales than the 
cold chemistry timescales, explaining the lack of cold chemistry tracers. Efficient photodesorption 
of CO into the gas phase would also explain the lack of N2H + and H2CO, while its impact on 
the deuterium fractionation is harder to assess. The same processes are probably present in disks 
around T Tauri stars as well, but because their disks are overall colder there is still enough material 
protected from vertical mixing and photodesorption on long enough timescales for large amounts 
of N 2 H+, DCO+ and H 2 CO to form. 

The lack of CH3OH detections does not put strong constraints on the CH3OH/H2CO abun- 
dance ratio, since H2CO is barely detected and the CH3OH transitions in this spectral region are 
more than an order of magnitude weaker than the observed H2CO transitions. To put stronger con- 
straints on CH3OH abundances in disks instead requires targeted observations of the most intense 
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CH.sOH lines. 



6. Conclusions 

Protoplanetary disks exhibit a rich chemistry that varies significantly between different objects 
within the same star forming region. Some of this variation can be understood in terms of the 
central star and its heating of the disk - the cold chemistry tracers N2H + , DCO + , DCN and 
H2CO are only detected toward T Tauri stars in our disk sample of four T Tauri stars and two 
Herbig Ae stars. Tracers of photochemistry, especially CN and HCN, show no clear dependence 
on quiescent stellar luminosity within the sample. Deuterium fractionation also seems to depend 
on parameters other than the disk temperature structure. For these chemical systems, the impact 
of other sources of irradiation, e.g. accretion shocks and X-rays, as well as the disk structure and 
grain characteristics may all be more important for the chemical evolution than the quiescent stellar 
luminosity. Investigating the relative importance of these different disk and star characteristics 
requires a combination of detailed modeling of the current sample, an increase in the number of 
sources to boost the statistics and span more parameters - especially a larger range of accretion 
rates and disks around the intermediate F stars - and targeted observations of rare isotopes of CO 
and HCO + to extract accurate abundance ratios. While the chemical evolution in protoplanetary 
disks is clearly complex, the qualitative agreement between at least parts of the early DISCS results 
and our current chemical understanding is promising for the ongoing modeling of these objects. 
The key results so far are listed below. 

1. Six disks in Taurus (DM Tau, AA Tau, LkCa 15, GM Aur, CQ Tau and MWC 480) have 
been surveyed for 10 molecules, CO, HCO+, DCO+, CN, HCN, DCN, H 2 CO, N 2 H+, CH 3 OH 
and C-C3H2, with a high detection rate and large chemical variability. 

2. The brightest molecular lines, CO 2-1, HCO+ 3-2, CN 3-2 and HCN 3-2 are detected toward 
all disks, except for HCN toward CQ Tau. Other molecular lines tracing different types of 
cold chemistry, N2H + , DCO + , DCN and H2CO, are only detected toward disks around T 
Tauri stars, indicative of a lack of cold regions around Herbig Ae stars for long enough time 
scales. 

3. Both the absolute CN flux and the CN/HCN ratio vary significantly among the observed 
disks and their variation seems independent of stellar luminosity, suggestive of that other 
parameters such as accretion luminosity and dust growth and dust settling play an important 
role for the chemical evolution in disks. 

4. Among the cold chemistry tracers the DCO + /N2H + ratio varies by an order of magnitude 
suggesting that the deuterium fractionation depends on other parameters, including the ra- 
diation field, beyond the amount of cold material present in the disk. 
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Fig. 1. — Spatially integrated spectra from DM Tau, AA Tau, LkCa 15 GM Aur, CQ Tau and 
MWC 480, displaying the chemical variation between the different disks. The line profiles of CO, 
HCO + , HCN and one of the CN lines have been scaled down for visibility. 
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Fig. 2. — Velocity channel maps for all detected molecules toward DM Tau . The first contours as well as 
the contour step values are ~ 2-cr, or 0.15 Jy beam -1 for CO and HCO + , 0.12 Jy beam.- 1 for N2H + and 
HCN, 0.07 Jy beam" 1 for DCO+, H 2 CO 4-3 and the CN lines, and 0.04 Jy beam- 1 for H 2 CO 3-2. he 
cross indicates the continuum (stellar) position. Note the different velocity resolution for the bottom four 
transitions. 



24 



1 




* 




0.1 




\ 




0.01 




•» 




0.001 


DM 


Tau 


♦ 



i 

0.1 


t 




0.01 




. 


0.001 


AA 


Tau 





* . 






• 




0.1 


•» 




0.01 


« 




0.001 


LkCe 15 


t 



1 
0.1 


h 






0.01 






•» 




0.001 


GM Aur 






* 



1 






» 


0.1 








0.01 






_ 


0.001 


co. 


Tau 









♦ 




1 
0.1 




\ 




0.01 








0.001 


MWC 


480 





0.1 1 10 

X [mm] 



DM Tau 
Continuum 



AA Tau 
Continuum 



LkCa 15 
Continuum 



< 5 



GM Aur 
Continuum 



CQ Tau 

Continuum 



sr\ 



Kc 480 
Continuum 




AA Tau 
CO 2-1 



85 deg 
4.5"/3,9" 




CQ Tau 
CO 2-1 



(S) 



-74 deg 
3.3"/2,7" 




DM Tau' : 
5\ HCO* 3-2 




AA Tau 
HCO* 3-2 



LkCa 15 ,'" 
HCO* 3-2 



GM Aur 
HCO* 3-2 



CQ Tau 
HCO* 3-2 



MWC 480 
5f HCO* 3-2 



DM Tau 

HCN 



5 0-5 5 

5 




5 0-5 

Aa ["1 



AA Tau 
HCN 



LkCa 15 

5[ HCN .;.-' 



5 0-5 5 

5 



GM Aur 

HCN 



^ 



5 0-5 a 

5 



CQ Tau 
HCN 



MWC 480 

HCN 








-5 5 

5 




AA Tau 
CN 



-5 5 

5 




-5 5 

5 



GM Aur 
CN 



O 




MWC 480 
CN 



Fig. 3. — Literature spectral energy distributions with our continuum points superimposed in red 
(col 1), 267 GHz continuum maps (col 2) and zero- and first moment maps of CO 2-1, HCO + 3-2, 
HCN 3-2 and CN 2-1 (col 3-6). In the moment maps the contours show the integrated intensity in 
steps of 0.8, 0.5, 0.32 and 0.32 Jy beam" 1 km s" 1 for CO, HCO+, HCN and CN respectively. The 
velocity gradient is defined over the same range for all molecules toward each source. The derived 
position angles and CO FHWM major and minor axes derived from fitting ellipses to the CO maps 
are shown in the lower right corner of each CO panel. 
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Fig. 4. — Integrated intensity line ratios and upper and lower limits for CN/HCN and DCO + /N2H + 
and H 2 CO/N 2 H + . 



Source 



Table 1. Central star and accretion data. 



RA 



DEC 



Sp. type L„ (L ) M (10" 9 M Q yr" 1 ) 



DM Tau 04 33 48.73 


+ 18 10 10.0 


M1.5 a 


0.3 a 


0.8 b 


AA Tau 04 34 55.42 


+24 28 53.2 


K7 C 


0.59 c 


7 C 


LkCa 15 04 39 17.78 


+22 21 03.5 


K3 a 


1.2 a 


8 b 


GM Aur 04 55 10.98 


+30 21 59.4 


K5 a 


0.93 a 


3 b 


CQ Tau 05 35 58.47 


+24 44 54.1 


F2 d 


8.0 C 


<8 f 


MWC 480 04 58 46.27 


+29 50 37.0 


A3 g 


11. 5 g 




a Espaillat et al. (submitt 


3d) 








iBersin et al.l J2004) 




4White & Ghed J200lh 










iThe et al.l jl994h 










^Mannings & Sargentl dl997l^ 









iGarcia Lopez et al.1 J2006I ') 
^Simon et all J2000h 
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Table 2. Disk data. 



Source PA co (deg) i co (deg) R CO (AU) R hoie (AU) M D (M Q ) 



DMTau 


-25 a 


-35 a 


890 a 


3 C 


0.02 d 


1.0° 


,1.5' 


AA Tau 


-89 s 


75 h 


995 g 




0.01 d 


0.3 C 


, 1.0' 


GM Aur 


46 b 


49 b 


445 b 


24 c 


0.03 d 


1.0° 


, 1.2' 


LkCa 15 


61 a 


52 a 


905 a 


58 1 


0.05 d 






CQ Tau 


-43J 


2f)J 


200 J 




0.009 k 


0.7J 


, 0.6 1 


MWC 480 


-32 a 


38 a 


740 a 




0.11 m 


0.8 n 


,1.4' 



iPietu et alj 1120071 ) - the reported PAs have been converted to conform with its conventional 
definition as the angle of the projected disk major axis. 

iDutrev et"aj] (2008) 

d Calvet et al.1 1120051 1 

^Andrews fc Williams! J2005h 

JRicci et af] ll2009h 

lAndrews fe Williams] J2007t) 

teessler-Silaccil J2004> 

^Kessler-Silaccil IJ2004I ') 

'Espaillat et al. (submitted) 

IChapillon et al.1 I I2008I ) — the reported PA has been converted to conform with its conventional 
definition as the angle of the projected disk major axis. 

^Mannings fc Sargent] feOOCh 

iTesti et all (1200311 

" IChiang et al.1 d200l l l 

■ JHamidouche et al.l J2OO6) 



Table 3. Spectral setup: 1.4 mm setting (214.670-218.653 GHz, LSB and 226.671-230.654 GHz, USB). 

Chunk Frequency range Channels Resolution Lines Frequency T up 

(GHz) (km s- 1 ) (GHz) (K) 

LSB 



S31 


216.064-216.168 


256 


0.56 


DCO+ 3-2 


216.1126 


21 


S18 


217.155-217.259 


256 


0.56 


DCN 3-2 


217.2386 


21 


S10 


217.811-217.915 


256 


0.56 


C-C3H2 61 6 — 5o 5 


217.8221 


39 


S06 


218.140-218.244 


61 


2.23 


H 2 C0 3 3-2o2 


218.2222 


21 



S03 


218.392-218.496 


256 


0.56 


CH3OH-E 4 2 1 - 3i 2 


218.4401 


46 


USB 


S01 


226.671-226.775 


128 


1.07 


CN2 2 2-ln 


226.6794 


16 


S03 


226.828-226.932 


256 


0.54 


CN 2 3 3/4/2 — I2 2/3/1 


226.8747 


16 


S47 


230.468-230.572 


256 a 


0.53 


CO 2-1 


230.538 


17 



toward LkCa 15 
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Table 4. Spectral setup: 1.1 mm setting (265.777-269.760 GHz, LSB and 277.778-281.761 GHz, USB). 



Chunk 


Frequency range 
(GHz) 


Channels 


Resolution 

(km s- 1 ) 


Lines 


Frequency 
(GHz) 


T 
- 1 - up 

(K) 


LSB 


S47 
S27 


265.858-265.962 
267.499-267.603 


256 
256 


0.46 
0.46 


HCN 3-2 
HCO+ 3-2 


265.8862 
267.5575 


26 
26 


USB 


S22 
S46 


279.499-279.603 
281.500-281.604 


256 
128 a 


0.44 
0.87 


N 2 H+ 3-2 
H 2 C0 4i4-3i3 


279.5117 
281.5269 


27 
16 



toward LkCa 15 



Table 5. SMA Observation Log. 



Date 


Sources 


Setting 3 - 


Baselines b (Antennas) 


7~225GHz 


Tsys(Kr 


2009 Nov 20 


LkCa 15, MWC 480 


1.1 mm 


8-64(7) 


0.04-0.08 


107-528 


Nov 22 


LkCa 15 


1.4 mm 


9-52(6) 


0.04-0.1 


74-103 


Nov 28 


MWC 480 


1.4 mm 


10-108(6) 


0.04-0.09 


107-528 


Dec 03 


DM Tau 


1.4 mm 


9-52(7) 


0.05-0.1 


87-125 


Dec 08 


DM Tau, CM Aur 


1.1 mm 


8-64(7) 


0.05-0.09 


116-316 


Dec 09 


GM Aur 


1.4 mm 


7-52(7) 


0.06-0.1 


87-216 


Dec 15 


AA Tau, CQ Tau 


1.1 mm 


9-64(6) 


0.06-0.09 


111-287 


Dec 19 


CQ Tau 


1.4 mm 


9-58(8) 


0.06-0.15 


79-195 


Dec 23 


AA Tau 


1.4 mm 


9-52(7) 


0.05-0.09 


73-150 


2010 Mar 22 


LkCa 15 


1.1 mm 


4-34(7) 


0.06-0.1 


83-151 



a See Tables [3] and [4] 

b Projected baseline in meters 

c Double sideband (DSB) system temperature 



Table 6. Beam sizes and position angles for each source and setting. 



Sources Setting Beam size (arc sec) Beam PA (degrees) 



DM Tau 


1.1 


mm 


3.3 X 2.6 


DM Tau 


1.4 


mm 


3.9 x 3.3 


AA Tau 


1.1 


mm 


2.9 x 2.5 


AA Tau 


1.4 


mm 


3.9 x 3.2 


LkCa 15 


1.1 


mm 


3.3 x 2.7 


LkCa 15 


1.4 


mm 


3.8 x 3.2 


GM Aur 


1.1 


mm 


3.2 x 2.7 


GM Aur 


1.4 


mm 


3.8 x 3.2 


CQ Tau 


1.1 


mm 


2.9 x 2.5 


CQ Tau 


1.4 


mm 


3.5 x 3.0 


1WC 480 


1.1 


mm 


3.2 x 2.5 


IWC 480 


1.4 


mm 


3.8 x 2.9 



19 

19 
27 
10 
26 
14 
26 
14 
27 
-70 
21 



28 



Table 7. Integrated spectra for DM Tau. Spectra for all sources are available on-line. 



CO 2-1 


HCO+ 3-2 DCO+ 3-2 N 2 H+ 3-2 H 2 CO 4-3 H 2 CO 3-2 HCN 3-2 


DCN 3-2 


CN 2 3 -l 2 CN 2 2 -li 


Vel. Flux 


Vel. Flux Vcl. Flux Vcl. Flux Vel. Flux Vcl. Flux Vel. Flux 


Vcl. Flux 


Vel. Flux Vel. Flux 


km/s Jy 


km/s Jy km/s Jy km/s Jy km/s Jy km/s Jy km/s Jy 


km/s Jy 


km/s Jy km/s Jy 



-8.50 

-8.00 

-7.50 

-7.00 

-6.40 

-5.90 

-5.40 

-4.80 

-4.30 

-3.80 

-3.30 

-2.70 

-2.20 

-1.70 

-1.10 

-0.60 

-0.10 

0.40 

1.00 

1.50 

2.00 

2.60 

3.10 

3.60 

4.10 

4.70 

5.20 

5.70 

6.30 

6.80 

7.30 

7.80 

8.40 

8.90 

9.40 

10.00 

10.50 

11.00 

11.50 

12.10 

12.60 

13.10 

13.70 

14.20 



0.09 

-0.07 

0.05 

0.19 

-0.23 

0.13 

-0.12 

-0.13 

-0.06 

0.01 

0.08 

0.08 

-0.06 

-0.30 

0.22 

-0.33 

0.31 

-0.12 

0.22 

-0.10 

-0.26 

0.07 

0.34 

0.29 

0.67 

0.60 

3.99 

7.73 

8.10 

4.84 

0.85 

0.18 

0.05 

0.49 

-0.13 

-0.07 

0.07 

0.28 

-0.11 

-0.29 

-0.04 

0.06 

0.26 

-0.05 



-8.70 

-8.20 

-7.80 

-7.30 

-6.90 

-6.40 

-6.00 

-5.50 

-5.10 

-4.60 

-4.20 

-3.70 

-3.20 

-2.80 

-2.30 

-1.90 

-1.40 

-1.00 

-0.50 

-0.10 

0.40 

0.90 

1.30 

1.80 

2.20 

2.70 

3.10 

3.60 

4.00 

4.50 

5.00 

5.40 

5.90 

6.30 

6.80 

7.20 

7.70 

8.10 

8.60 

9.00 

9.50 

10.00 

10.40 

10.90 



0.03 

0.20 

0.09 

-0.63 

0.01 

-0.09 

0.22 

0.35 

0.06 

-0.39 

-0.19 

0.49 

-0.06 

0.04 

-0.49 

-0.25 

-0.08 

0.03 

-0.09 

-0.08 

0.48 

-0.07 

-0.36 

-0.01 

0.44 

0.28 

-0.19 

0.13 

-0.12 

0.19 

0.96 

2.32 

2.34 

3.11 

2.36 

0.81 

-0.01 

0.41 

-0.30 

-0.32 

0.15 

0.53 

0.15 

-0.22 



-8.90 

-7.80 

-6.70 

-5.50 

-4.40 

-3.30 

-2.20 

-1.00 

0.10 

1.20 

2.30 

3.50 

4.60 

5.70 

6.90 

8.00 

9.10 

10.20 

11.40 

12.50 

13.60 

14.70 

15.90 

17.00 

18.10 



-0.06 

0.07 

-0.02 

-0.03 

-0.02 

0.08 

0.04 

-0.02 

-0.02 

-0.14 

0.06 

0.04 

0.10 

0.26 

0.24 

-0.01 

0.05 

0.02 

0.05 

-0.02 

0.04 

-0.05 

-0.07 

0.00 

-0.03 



-8.70 

-7.80 

-7.00 

-6.10 

-5.20 

-4.30 

-3.50 

-2.60 

-1.70 

-0.90 

0.00 

0.90 

1.80 

2.60 

3.50 

4.40 

5.20 

6.10 

7.00 

7.90 

8.70 

9.60 

10.50 

11.30 

12.20 

13.10 

14.00 

14.80 

15.70 

16.60 

17.40 

18.30 



0.07 

-0.06 

0.09 

0.09 

-0.31 

0.12 

-0.19 

0.02 

0.08 

-0.01 

0.31 

-0.02 

-0.21 

0.05 

0.07 

-0.02 

0.47 

0.31 

0.34 

-0.03 

0.05 

-0.03 

0.20 

-0.09 

-0.01 

-0.00 

-0.01 

0.18 

-0.05 

-0.06 

-0.05 

0.06 



-8.80 

-7.90 

-7.00 

-6.20 

-5.30 

-4.40 

-3.60 

-2.70 

-1.80 

-1.00 

-0.10 

0.80 

1.60 

2.50 

3.40 

4.20 

5.10 

5.90 

6.80 

7.70 

8.50 

9.40 

10.30 

11.10 

12.00 

12.90 

13.70 

14.60 

15.50 

16.30 

17.20 

18.10 

18.90 



-0.01 

0.03 

0.01 

-0.20 

-0.10 

0.23 

0.07 

0.06 

0.02 

0.24 

0.04 

0.08 

-0.10 

0.01 

-0.22 

-0.17 

0.01 

0.46 

0.09 

-0.02 

0.03 

-0.02 

-0.12 

-0.04 

-0.04 

-0.20 

-0.16 

0.07 

-0.09 

0.12 

0.07 

-0.10 

-0.01 



-7.70 

-5.50 

-3.20 

-1.00 

1.20 

3.40 

5.70 

7.90 

10.10 

12.40 

14.60 

16.80 



-0.00 

0.01 

0.03 

0.03 

0.00 

0.03 

0.10 

0.02 

-0.03 

-0.00 

0.05 

-0.06 



-8.70 

-8.30 

-7.80 

-7.40 

-6.90 

-6.40 

-6.00 

-5.50 

-5.10 

-4.60 

-4.20 

-3.70 

-3.20 

-2.80 

-2.30 

-1.90 

-1.40 

-0.90 

-0.50 

0.00 

0.40 

0.90 

1.30 

1.80 

2.30 

2.70 

3.20 

3.60 

4.10 

4.60 

5.00 

5.50 

5.90 

6.40 

6.80 

7.30 

7.80 

8.20 

8.70 

9.10 

9.60 

10.00 

10.50 

11.00 



-0.04 

-0.01 

-0.14 

0.16 

0.04 

-0.01 

0.58 

-0.55 

0.21 

0.06 

0.18 

0.16 

-0.31 

-0.17 

-0.31 

-0.29 

0.03 

0.10 

-0.31 

-0.07 

0.45 

-0.27 

-0.06 

-0.07 

0.17 

-0.21 

0.02 

-0.04 

0.03 

0.13 

0.82 

1.20 

1.29 

1.57 

0.82 

0.37 

0.46 

-0.17 

-0.14 

-0.05 

-0.16 

0.18 

0.33 

0.20 



-8.50 

-7.40 

-6.20 

-5.10 

-4.00 

-2.90 

-1.80 

-0.60 

0.50 

1.60 

2.70 

3.90 

5.00 

6.10 

7.20 

8.30 

9.50 

10.60 

11.70 

12.80 

13.90 

15.10 

16.20 

17.30 

18.40 



-0.04 

-0.06 

0.02 

0.03 

-0.01 

-0.04 

-0.04 

0.04 

0.08 

0.06 

-0.01 

-0.10 

-0.02 

0.09 

-0.11 

0.12 

-0.05 

0.05 

-0.01 

-0.06 

0.01 

-0.02 

-0.06 

0.00 

-0.05 



-8.60 

-8.10 

-7.50 

-7.00 

-6.40 

-5.90 

-5.40 

-4.80 

-4.30 

-3.80 

-3.20 

-2.70 

-2.20 

-1.60 

-1.10 

-0.50 

0.00 

0.50 

1.10 

1.60 

2.10 

2.70 

3.20 

3.80 

4.30 

4.80 

5.40 

5.90 

6.40 

7.00 

7.50 

8.00 

8.60 

9.10 

9.70 

10.20 

10.70 

11.30 

11.80 

12.30 

12.90 

13.40 

14.00 

14.50 



-0.00 

0.04 

0.09 

0.06 

-0.01 

0.29 

-0.00 

-0.13 

0.04 

0.11 

0.05 

-0.15 

0.12 

0.04 

-0.12 

-0.05 

-0.08 

-0.06 

-0.11 

-0.19 

0.14 

0.09 

-0.04 

0.16 

0.51 

0.90 

1.59 

1.45 

1.67 

1.73 

0.65 

0.14 

0.06 

0.02 

0.08 

0.05 

-0.16 

0.03 

0.02 

0.22 

-0.01 

-0.01 

0.11 

-0.07 



-8.60 

-7.50 

-6.40 

-5.30 

-4.30 

-3.20 

-2.10 

-1.10 

0.00 

1.10 

2.20 

3.20 

4.30 

5.40 

6.50 

7.50 

8.60 

9.70 

10.80 

11.80 

12.90 

14.00 

15.10 

16.10 

17.20 



-0.02 

-0.12 

-0.07 

-0.17 

0.01 

-0.05 

-0.18 

-0.00 

-0.02 

0.14 

-0.05 

0.12 



20 
12 
22 
00 
0.08 
03 
07 
01 
06 
00 
08 
06 
08 
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Table 7- 


— Continued 
















CO 2-1 


HCO+ 3-2 DCO+ 


3-2 N 2 H+ 3-2 H 2 CO 4-3 H 2 CO 3-2 


HCN 3-2 


DCN 3-2 


CN 2 3 -l 2 CN 2 2 -li 


Vel. Flux 


Vcl. Flux Vcl. Flux Vcl. Flux Vcl. Flux Vcl. Flux 


Vcl. Flux Vel. Flux 


Vel. Flux Vcl. Flux 


km/s Jy 


km/s Jy km/s Jy km/s Jy km/s Jy km/s Jy 


km/s Jy km/s Jy 


km/s Jy km/s Jy 


14.70 -0.20 


11.30 -0.22 . 


















11.40 -0.03 . 






15.00 0.05 . 






15.20 -0.25 


11.80 0.01 . 


















11.90 -0.01 . 






15.60 0.24 . 






15.80 0.09 


12.20 0.13 . 


















12.30 -0.12 . 






16.10 0.18 . 






16.30 0.14 


12.70 0.04 . 


















12.80 0.25 . 






16.60 -0.13 . 






16.80 0.11 


13.10 -0.14 . 


















13.30 0.23 . 






17.20 -0.31 . 






17.40 0.44 


13.60 -0.36 . 


















13.70 0.01 . 






17.70 -0.18 . 






17.90 -0.08 


14.10 -0.05 . 


















14.20 -0.16 . 






18.20 -0.02 . 






18.40 0.37 


14.50 0.43 . 


















14.60 -0.14 . 






18.80 0.02 . 






18.90 -0.38 


15.00 -0.11 . 


















15.10 0.31 . 


















15.40 0.23 . 


















15.50 0.52 . 


















15.90 0.30 . 


















16.00 -0.34 . 


















16.30 -0.31 . 


















16.50 0.00 . 


















16.80 0.26 . 


















16.90 0.42 . 


















17.20 -0.43 . 


















17.40 0.34 . 


















17.70 -0.05 . 


















17.80 0.35 . 


















18.20 -0.22 . 


















18.30 0.33 . 


















18.60 0.00 . 


















18.80 0.02 . 















Table 8. Continuum flux in Jy and spectrally integrated line intensities and 1-a upper limits in 

Jy km s _1 integrated over the disk area. 



Species 



DM Tau 



AA Tau 



LkCa 15 



GM Aur 



CQ Tau 



MWC 480 



218 GHz cont. 
267 GHz cont. 

CO 2-1 
HCO+ 3-2 
DCO+ 3-2 
N 2 H+ 3-2 
H 2 CO 4-3 
H 2 CO 3-2 
HCN 3-2 
DCN 3-2 
CN 2 3 -l 2 
CN 2 2 -li 



0.104[0.010] 
0.143[0.014] 

14.87[0.12] a 
5.34[0.12] a 
0.71[0.04] 6 
0.97[0.07] c 

(0.29[0.11]) c 
0.35[0.03] c 
2.92[0.11] 6 

<0.10 a 
4.74[0.08] a 
0.57[0.07]° 



0.073 [0.007] 
0.113[0.011] 

8.20[0.16] a 
2.32[0.19] 6 

<0.18 a 
(0.66[0.16]) c 

<0.52 i> 
(0.16[0.08]) 6 
0.59[0.13] c 

<0.18 a 

1.76[0.09] b 

[0.08] b 



0.130[0.013] 
0.214[0.021] 

13.94[0.15] a 
5.19[0.17] a 
0.51[0.08] 6 
0.71[0.14] c 
1.12[0.22] 6 
0.66[0.06] b 
5.52[0.21] a 
0.41[0.05]<> 
6.80[0.11] a 
1.16[0.12]° 



0.177[0.018] 
0.251[0.025] 

19.41[0.11] a 
5.01[0.17] a 

<0.15 a 
1.37[0.14] a 
0.57[0.23] 6 
0.56[0.08] 6 
1.27[0.12] c 

<0.09 a 
0.96[0.08] a 

<0.19 a 



0.150[0.015] 
0.258[0.026] 

3.10[0.18] a 
0.75[0.24] a 
<0.23 a 
<0.35 a 
<0.42 a 
<0.21 a 
<0.50 a 
<0.17 a 
0.22[0.11] b 
<0.22 b 



0.255[0.026] 
0.427[0.043] 

22.00[0.20] a 
4.65[0.31] a 
<0.31 a 
<0.75 a 
<0.78 a 
<0.28 a 
2.34[0.19] 6 

<0.22 a 
3.29[0.14] a 
0.61[0.15] a 



The integrated intensities are extracted from a the complete CO disk, b 
uncertainties in the integrated intensities are in square parentheses. 



3/4 of the CO disk and c 1/2 of the CO disk. The 1- 
Intcnsities in parentheses are tentative detections. 



